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ABSTRACT 

We report a centimetre-wave (cm-wave, 5-31 GHz) excess over free-free emission 
in PNe. Accurate 31 and 250 GHz measurements show that the 31 GHz flux densities in 
our sample are systematically higher than the level of optically thin free-free continuum 
extrapolated from 250 GHz. The 31 GHz excess is observed, within one standard 
deviation, in all 18 PNe with reliable 31 and 250 GHz data, and is significant in 9 
PNe. The only exception is the peculiar object M 2-9, whose radio spectrum is that 
of an optically thick stellar wind. On average the fraction of non-free-free emission 
represents 51% of the total flux density at 31 GHz, with a scatter of 11%. The average 
31-250 GHz spectral index of our sample is (ofgf^) = — 0.43±0.03 (in flux density, with 
a scatter of 0.14). The 31-250 GHz drop is reminiscent of the anomalous foreground 
observed in the diffuse ISM by CMB anisotropy experiments. The 5-31 GHz spectral 
indices are consistent with both flat spectra and spinning dust emissivities, given the 
10% calibration uncertainty of the comparison 5 GHz data. But a detailed study of 
the objects with the largest cm-excess, including the low frequency data available in 
the literature, shows that present spinning dust models cannot alone explain the cm- 
excess in PNe. Although we have no deflnitive interpretation of our data, the least 
implausible explanation involves a synchrotron component absorbed by a cold nebular 
screen. We give flux densities for 37 objects at 31 GHz, and for 26 objects at 250 GHz. 

Key words: radiation mechanisms: general planetary nebulae: general radio contin- 
uum: ISM sub-millimetre 



1 INTRODUCTION 

An increasing amount of evidence supports the existence of 
a new continuum emission mechanism in the diffuse inter- 
stellar medium (ISM) at 10-30 GHz, other than free-free, 
synchrotron, or an hypothetical Raylcigh- Jeans tail of cold 



dust grains (e.g., Leitch et al.|l997, |Finkbeiner|2004 Wat- 
|son et aT]|2005[ ). As proposed by |Draine" fc Lazarianf ( |1998| ) 

a promising candidate emission mechanism is electric dipole 
radiation from spinning very small grains (VSGs), or spin- 
ning dust. Cosmic Background Imager (CBI) observations of 
the dark cloud LDN 1622 linked, on morphological grounds, 
the cm- wave emitters to the VSGs ( Casassus et al.|2006 l. 

The SED of the Helix nebula (NGC 7293, an evolved 
planetary nebula, PN) also exhibits a 31 GHz excess over the 
free-free continuum expected from 250 GHz measurements 
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( [Casassus et al.|2004[ ). Resolved 31 GHz images of the Helix, 
obtained with the CBI, provided morphological evidence for 
a new radio continuum component. 

PNe have traditionally served as test laboratories for 
nebular astrophysics. Their spectra are the product of the 
photoionisation by a single exciting star of its surrounding 
point-symmetric nebula, in a characteristic radial and ho- 
mologous expansion. PNe are archetypical free-free sources; 
the interpretation of their < 14 GHz continuum in terms of 



free-free emission is paradigmatic (e.g. Siodmiak & Tylenda 
|2001[ [Condon fc Kaplan|1998[ ) . Yet inconsistencies with the 
free-free paradigm exist in the literature, although not recog- 
nised as such. Hoare et al. (19921, in an investigation of 
the cm- to mm-wave SED of PNe, ascribed deviations from 
free-free spectra to instrumental uncertainties. [Garay et al.[ 
( 1989 1 compared continuum and radio-recombination line 



maps, finding systematic differences interpreted as non-LTE 
amplification or large Te variations. 
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Here we report on a survey of the cm- and mm-wave 
continuum in PNe acquired with the Cosmic Background 



Table 1. Observed and free-free spectral indices at Te = 7000 K. 



Imager (CBI, |Padin et al.||2002|) and with the Ses t IMag- 
ing Bolometer Array (SIMBA, |Nyman et aLl|2001[ ), at the 
Swedish-ESO Submillimetre Telescope (SEST). Th e FN 
sample is mostly selected from Casassus et al.| ( |2001 1, and 
represents the population of compact and mid-IR-bright 
PNe (it is not a complete sample). Section |2] discusses FN 
electron temperatures and relevant free-free spectral indices. 
Section [3] and [4]present the CBI and SIMBA data, respec- 
tively. Section [5] compares the two datasets, which are in 
turn compared to previous measurements in Section [6] Sec- 
tion [7] interprets our findings in terms of possible emission 
mechanisms. Section|8]analyses the SEDs of individual FNe. 
Section[9]concludes. All the data analysis was carried out us- 
ing Ferl Data Language (http://www.pdl.org), unless oth- 
erwise stated. The error bars on all Figures are one standard 
deviation (±1 a, the error bars have a total length of 2 a). 



2 ADOPTED PLANETARY NEBULA 
ELECTRON TEMPERATURES AND 
FREE-FREE SPECTRAL INDICES 

The radio-frequency free-free spectral index is roughly con- 
stant at a ~ —0.15 in flux density, F{i') — F{vo){p / i'o)°' ■ 
Yet a has a residual dependence on frequency, electron tem- 
perature To, and Helium contribution. Since our aim is to 
perform a precision analysis of the FN SEDs, we calculated 
accurate free-free indices for the frequencies relevant to this 
work. The values we obtained are exact for non-relativistic 
electrons (i.e. Te < 10* K). Full details are given in Sec, [a] 
The electron temperature in FNe has been the matter 



of significant debate. Feimbert ( 1971 1 first observed that 



optical recombination lines (ORLs) and coUisionally-excited 
lines (CELs) give broadly different values. Since then, sev- 
eral new Te diagnostics have been reported, all pointing to- 
wards the conclusion that ORLs and continuum diagnostics 
give systematically lower values than [O ill] A4363/A4959, 
by a factor of ~ 2 (and up to 10). 



The current view, as summarised by Liu ( 2006[ ), is that 
of a bi-abundance model in which FNe (and to a lesser extent 
Hll regions) contain dense, hydrogen-deficient (with metal- 
icity Z ~ 1) and cool inclusions accounting for most of the 
ORL emission. In the bi-abundance model the bulk of the 
nebula, by mass, is best represented by the CEL diagnostics, 
with Te ~ 10000 - 15000 K. In the case of NGC 6153, the 
free-free luminosity of the ORL emitting gas is ~1/10 that 
of the CEL emitting gas, and the hydrogen mass in the ORL 
gas is only ~ 1 — 2 % that of the CEL gas (Mike Barlow, 
private communication, and |Liu et aLl|2000| their model 
IH3). Assuming the ORL Te is ~ 1/2 the CEL Te, an aver- 
age Te, weighted by free-free luminosity gives, in the case of 
NGC 6153, {Te) = 0.95re(CEL). 

Our strategy was to adopt a single conservative value 
for the whole FN sample. We thus adopted Te = 7000 K, a 
value that approaches a canonical Te for Hll regions (e.g., 
|Dickinson et al. 2003) , but is unrealistically cold for FNe. 
This choice minimises free- free deviations in the dataset, as 
cooler Te give steeper 31-250 GHz spectral indices. Since the 
emissivity of the He"*" continuum gives somewhat steeper 
indices than for H, we also included a contribution from 





„36 


„250 
"31 




Free-free indices 




-0.125'" 


-0.137 


-0.152 




PN sample 




-0.06 ± 0.01(0. 22)'''= 


-0.18 ±0.03(0.21)'* 


-0.43 ±0.03(0.14)"= 



" Tlie free-free indices bear a Icr (one standard deviation) 
uncertainty of 10~^, stemming from the uncertain contribution 
of the He"*" free-free continuum. 

The root-mean-square (rms) scatter is indicated in parenthesis. 

From all 31 objects with 5 GHz flux densities and X^gj < 2. 

From all 32 objects with Xcbi ^ ^' 
" From all 17 objects with 250 GHz flux densities and Xcbi ^' 



the He^ continuum, assuming a large but conservative 10% 
abundance of He^""" relative to H"*". Table [T] summarises the 



free- free spectral indices relevant to our analysis (see Sec. A3 
for details). 



3 CBI OBSERVATIONS 

The CBI is a planar interferometer array with 13 antennas, 
each 0.9 m in diameter, mounted on a 6 m tracking plat- 
form, which rotates in parallactic angle to provide uniform 
iiw-coverage. Its uniform-weight synthesised beam is ~6 ar- 
cmin. The CBI receivers operate in 10 frequency channels, 
with 1 GHz bandwidth each, giving a total bandwidth of 
26-36 GHz. It is located in Llano de Chajnantor, Atacama, 
Chile. In this Section we describe the analysis of total inten- 
sity and polarisation observations of FNe. 



3.1 Total intensity 

We acquired CBI total intensity observations of a sample of 
36 FNe, whose flux densities are hsted in Table [2] The visi- 
bility data were reduced and edited using a special-purpose 
package (CBICAL, developed by T.J. Pearson). Flux cal- 
ibration was obtained by reference to Mars or Jupiter, 
and tied to the WMAP temperature for Jupiter, Tjup = 
147. 3± 1.8 K ( Readhead et al._2004^ and references therein). 
The fractional uncertainty on Tjup, o"(Tjup) = 1.2 %, affects 
the overall calibration of the CBI. 

Approximate cancellation of ground and Moon contam- 
ination was obtained by differencing with a reference field 
at the same declination but offset in hour angle by the du- 
ration of the on-source integration. In general we used an 
on-source integration time of 8 min. We found that for point 
sources the acquisition of a reference field is not mandatory, 
since restricting to the longer baselines avoids any ground 
contamination (e.g., Casassus et al. 20061. We thus used 



the differenced dataset when available, and discarded base- 
lines shorter than 200 A otherwise. The baseline-restricted 
and differenced datasets have similar signal-to-noise ratios 
on point source measurements. 

We extracted integrated flux densities by fitting 
parametrised models to the visibilities. The uncertainties 
on the inferred flux densities correspond to the Ax^ = 1 
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contour. Point source models gave better reduced val- 
ues than elliptical gaussians or core-halo models, except for 
NGC 1360, which we fit with a 4.9 x 2.3 arcmin^ eUipti- 
cal Gaussian (full width half maximum, FWHM), with its 
major axis at a position angle of 22.9 deg, positive East of 
North. 

The use of point-source models to extract flux densities 
assumes there is no contaminating emission within the CBI 
primary beam, which has a FWHM of 45.2 arcmin. Addi- 
tional sources within the primary beam result in high re- 
duced values (xcbi hereafter). We checked by inspection 
that extended emission is present for Xcbi greater than ~2.0, 
which is the cutoff we adopted to select reliable point source 
flux densities. We attempted to extract flux densities from 
images of the flelds with extended emission, and obtained 
the same value as listed in Table [2] The extended emission 
fills the CBI primary beam, is featureless and difficult to 
separate from the PN. In addition to the extended emission, 
there are also bright point sources near NGC 6537. 19W32 
is altogether drowned by the diffuse background. Thus the 
sky-plane measurements have the same degree of unreliabil- 
ity as the visibility fits. 

A measurement is considered reliable if its uncertainty, 
as estimated from the point-source models, indeed corre- 
sponds to a 1 (T deviation. Several integrations on different 
nights for a subsample of 22 objects confirmed that, when 
Xcbi < 2.0, all nightly measurements are consistent within 
3 a with the fiux densities derived from model-fitting the 
combined visibility dataset. We also checked for consistency 
between the differenced and baseline-restricted datasets. 

In addition to the statistical uncertainties due to ther- 
mal noise, the final error bars also include uncertainties 
in the overall calibration of the CBI for the PNe obser- 
vations. The systematics stem in part from uncertainties 
in our primary calibrator (Jupiter), and in part from gain 
variations during a given night. We thus compared PN flux 
densities obtained with two different calibrators, during the 
same night. We measured the fractional flux density differ- 
ences. Again = 1-Fcaii --Fcai2l/-Fcaii, whcu Calibrating against 
calibrator 1 or 2. On the night of 15-Jan-2003, NGC 3918 
was observed along with Jupiter and Tau A; on 3-Jun-2001 
NGC3242 was observed along with J1230+123 and Mars. 
Using the same reduction (CBICAL) scripts and analysis 
tools as those used to produce the flux density reported in 
Table[2] we obtain Again = 0.68% for NGC 3918, and 0.80% 
for NGC 3242. For the PNe observations, the uncertainty on 
the CBI response is thus of order cTgain = 1%. 

Including the uncertainty on Tjup, we take as a con- 
servative estimate of fractional systematic uncertainties 
CTaystematic = 2%. The flual Uncertainty on the flux densities 
-FcBi reported in Table[2]is a^ot^i = (iTsystematic -Fcbi)^ + o"^, 
where a derives from model-fltting the combined visibility 
datasets. 

The model-flt spectral indices over 26-36 GHz are all 
consistent with free-free emission (see Table [l]). The ac- 
curacy of the CBI indices reaches ~Q.l for NGC 6369, 
NGC 6572, NGC 7009, IC 418, and NGC 3242. 

The CBI 31 GHz and comparison 5 GHz flux densities 
are also consistent with optically thin free-free, within the 
10% calibration uncertainty of the literature data (Table [T]). 
The only exceptions are Pe 1-7, M 2-9, SwSt 1, and Vy 2-2, 



Table 3. CBI 31 GHz polarisation measurements for four PNe. 
Upper limits are 99% C.L.. 



PNe 


Polarisation fraction 


NGC 7009 


< 8.5 % 


NGC 1360 


< 8.5 % 


NGC 246 


< 17 % 


NGC 7293 


< 2.4 % 



all with excesses at 31 GHz and point-source models with 
Xcbi < 2. 

In a comparison with the 14 GHz dataset, we obtain 
(qs'*) = 0.00 ±0.03 with an rms scatter of 0.19, and (all) = 
—0.23 ±0.03 with a scatter of 0.37. (afi) is steeper than the 
free-free value by 3.3 u, if the uncertainties on the 14 GHz 
are indeed smaller than 10% rms. 



3.2 CBI polarization upper limits. 

Four PNe were observed by the CBI in polarization dur- 
ing November 2002 and June 2005: NGC7009, NGC1360, 
NGC246 and NGC7293 (the Helix nebula). The polariza- 
tion capability of the CBI has been demonstrated with deep 
observations of CMB (Readhead et al. 20041 and also in 



Hi: regions ( [Cartwright et al.||2005| [Dickinson et al.|[2006 l 



Details of the polarization calibration and data reduction 
procedures can be found in the aforementioned papers. 

Natural-weighted maps of Stokes Q and U were made 
using the difmap package and comb ined to make polariza- 
tion intensity maps, P = + . All maps were consis- 
tent with noise - no signiflcant polarization was detected in 
any of the four PNe. Table |3] summarizes the results and 
includes the 3a upper limit on the polarization fraction, 
p = P/I, derived from the data. 

The 99% confldence level (CL) upper limit on the polar- 
ization fraction was calculated assuming the Q and U maps 
have normal distribution with zero mean and a standard 
deviation of a. The polarized intensity, P = \/ Q'^ + U'^ , 
follows the Rayleigh distribution f{P). The CL F{Pu) as- 
sociated to an upper limit P„ on the polarized intensity is 
thus F{P„) = /jj^" f{P)dP, and a CL of 99% corresponds to 
P„ « 3(7. Finally the upper limit on the fraction of polarized 
intensity is Pu/7max, where /max is the peak total intensity 
(which in units of Jy beam~^ is equal to the source flux 
density for unresolved sources). 

Strong radio polarisation is not expected for PNe 
since the emission is typically dominated by free-free emis- 
sion, which is intrinsically unpolarised. However, the excess 
anomalous component could perhaps be signiflcantly po- 
larised, depending on the physical mechanism producing the 
emission. The most constraining and interesting result here 
is for the Helix nebula where significant excess emission was 



observed at 31 GHz with the CBI (Casassus et al. 2004 1. 
Since free-free radiation is not polarised, we can limit the 
polarisation fraction on the excess emission in the Helix to 
<3.8-12 %, for 36-80 % free-free emission at 31 GHz. 
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Table 2. Summary of the CBI and SIMBA measurements. Reduced values, Xcbp ''^f^'' to parametrised models of the CBI visibilities. 
All flux densities are reported in mjy. 
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5 GHz and 14 GHz measurements from the compilation of [Acker et al.| | |l992[ |, taken mostly from the Parkes 64m data of |Milne &: | 
|Aller| | [l982| , and bearing a ^10% calibration uncertainty. 

"A is the difference between the observed 250 GHz flux density and optically thin free-free emission extrapolated from 31 GHz, if 
available, or 5 GHz otherwise. 

<5 = 100 X A/(i^(250 GHz) — A) is the percent fraction of 31 GHz emission not due to free-free emission. 

The 5 GHz flux for NGC 1360 is contaminated by a point source 
' Upper limits are 3 cr. 

4 SIMBA OBSERVATIONS 

SIMBA, at SESlQ is a 37-channel bolometer array, oper- 
ating at 1.2 mm (250 GHz). The half-power beam-width of 
a single element is 24". We observed each PN with a scan- 
ning speed of 80 " s~^, in 51 steps of 8 " orthogonal to 
the scan direction, thus obtaining 900" x 408" maps. The 
SIMBA scans are reduced with the standard recipe for point 



^ the SEST is situated on La Silla in Chile 



sources using the MOPSI package written by Robert Zylka 
(IRAM, Grenoble). Flux calibration is carried out by com- 
parison with Uranus maps. We extracted the calibration flux 
densities of Uranus by integrating the Uranus specific inten- 
sity map in a circular aperture 1 arcmin in radius about the 
centroid of an elliptical Gaussian fit. A residual sky back- 
ground was estimated from the median intensity in a ring 
1 arcmin wide and immediately surrounding the flux extrac- 
tion aperture. 

We observed 35 PNe during 3 different observing runs 
in 2001 and 2002. Flux densities were measured by integral- 
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Figure 1. SIMBA map of NGC 6302. The gray scale is in units of 
MJy/sr, and the peak intensity is 138 MJy/sr (or 2.1 Jy/beam, for 
SEST's circular beam with a FWHM of 24 arcsec at 250 GHz), x- 
axis and j/-axis show offset J2000 RA and DEC in arcsec. Contour 
levels at -4.9, -3.6, 0.2, 7.5, 19.5, 37.4, 62.4, 95.7, 138.4, MJy sr'^ . 
The SIMBA beam is indicated on the lower-left corner of the map. 



ing each PN specific intensity map in a circular aperture 
1 arcmin in radius about the centroid of an elliptical Gaus- 
sian fit to the PN image, and subtracting a residual sky 
background, as in the case of Uranus. 

The uncertainty assigned to each individual integration 
is the rms noise in a circular annulus surrounding our extrac- 
tion aperture (and taking into account correlated pixels). Up 
to three integrations from different days on individual ob- 
jects were combined to reduce the uncertainties. 

A ~10% calibration uncertainty affects SIMBA flux 
densities. The systematic error comes from the uncertainty 
in the brightness temperature estimate of Uranus and the 
uncertainty in the determination of telescope parameters, 
such as the aperture efficiency. A fractional systematic un- 
certainty with a standard deviation of 10% was therefore 
added in quadrature to the combined statistical uncertain- 
ties derived from the aperture photometry. 

Fig. [1] illustrates our SIMBA images. AU the SIMBA 
images can be approximated to elliptical Gaussians. The 
only case where the nebular solid angle Qn, as inferred 
from the Gaussian flts, is larger than 1.5 times the SIMBA 
beam, Jig, is NGC 6369, with fijv = 1.6r2s. In the case of 
NGC 6302 (shown on Fig. [TJ low level emission can be re- 
solved by SIMBA. The NE-SW low-level extension in Fig.[l] 
follows the bipolar axis of NGC 6302. The stripes and neg- 
atives in the SIMBA maps are proportional to the object 
flux densities. They are artifact inherent to the technique 
of fast scanning, which involves a high-pass fllter to cancel 
sky emission. These artifacts are also present on the Uranus 
calibration images, and are taken into account by the flux 
calibration procedure. 



5 CBI-SIMBA COMPARISON 

As summarised in Table [2] we have obtained 31 GHz and 
250 GHz flux densities for a total of 21 object, of which 3 
have XcBi > 2. Optically thin free-free emission at 31 GHz 



can be extrapolated to 250 GHz with a spectral index of 
a = -0.15. We find that 20 out of 21 objects have 250 GHz 
data consistent, within 1 a, with a deficit over the free-free 
level extrapolated from 31 GHz. Only 4 objects have < 1 a 
excesses at 250 GHz rather than a deficit (i.e. A > 0, in the 
notation of Table [2|. The 250 GHz deficit is significant in 9 
PNe, at 3 a. The only significant detection of an excess at 
250 GHz is in M 2-9, with a 5 cr excess at 250 GHz. 

The CBI-SIMBA spectral index, averaged over the 17 
measurements of aif listed in Table [2] that have Xcbi < 
2, and excluding M 2-9, is (0:3!°) = -0.43 ± 0.03, with a 
weighted rms scatter of 0.14. Thus, the observed (031°) is 
significantly different from the free-free value. 

We searched for a correlation between the 31 GHz 
excess and the mid-IR continuum due to VSGs, as ob- 
served in the anomalous CMB foreground. Unfortunately 
the IRAS 12/im band is contaminated by strong ionic lines 
as well as free-free emission, while ground-based lOfim spec- 
tra available from the literature are extracted from varying 
apertures. We can still compare the 31 GHz excess to the 
IRAS 25, 60, and 100 ^lm fiux densities. The IRAS 60^lm 
band turned out to give the best correlation with the 
250 GHz deficit. Fig. [2] plots the difference between the ex- 
pected level of free-free emission at 250 GHz, extrapolated 
from 31 GHz, against IRAS 60/im. The lO/xm dust emis- 
sion features compiled from the literature are indicated by 
'C' for PAHs, 'c' for SiC, 'O' for silicates, for weak and 
featureless continuum, and '*' when no data are available 
(as in Casassus et al.pOOl l. A straight line fit to the data, 
required to cross the origin, gives a dimensionless slope of 
(-3.05 ±0.18) 10"^. The cross-correlation with IRAS 60^m 
could simply reflect an overall scaling of the flux densities 
with nebular luminosity. The cm-excess does not seem re- 
lated to the 10 nm dust emission features. 

We also calculated the 250 GHz deflcit when extrapo- 
lating the free-free level from the 5 GHz data. The resulting 
deflcits are still significant, albeit for a reduced number of 
objects. The reduced signiflcance may be due in part to the 
10% calibration uncertainty of the 5 GHz data, but also to a 
5-31 GHz spectral index greater than the optically thin value 
(see Sec. |3|. In this case a straight line fit to the 250 GHz 
deficit as a function of IRAS 60//m band flux density, for the 
same objects as for the 31 GHz cross-correlation of Fig. [2] 
gives a dimensionless slope of (—2.65 ± 0.27) 10""^. 

Relative to the total 31 GHz flux density, the excess 
emission over the level of free-free emission extrapolated 
from 250 GHz represents ~0.3-0.7. In order to avoid in- 
duced correlations from scaling effects, we plot in Fig. [3] the 
ratio of the 31 GHz excess to the total 31 GHz flux density 
against the ratio of the IRAS 60 nm flux density over the 
5 GHz flux density. We have excluded from Fig. |3] all objects 
with large 1 a uncertainties on the fractional excess (> 0.3). 
The remaining 10 objects, not including M2-9, have an av- 
erage fractional excess of 0.51 ± 0.03 with an rms scatter of 
0.11. Interestingly, the fractional excess seems to correlate 
with the ratio of far-IR and radio flux, as expected for a 
dust-related emission mechanism at 31 GHz. 

We caution that the IR-radio slopes derived in this Sec- 
tion cannot be directly compared with the slopes observed 
in the ISM at large. PNe have much hotter dust tempera- 
tures ( 100-200 K, e.g., [Casassus et~al]|2001| ) than dark or 
cirrus clouds ( 15-25 K). 
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Figure 2. Deficit at 250 GHz from the expected level of free-free emission extrapolated from the CBI measurements, excluding M 2-9 
for clarity. The y— axis shows the difference between the SIMBA flux density and the free-free level extrapolated from 31 GHz, and the 
X-axis shows IRAS GOfim flux density. The dotted line is a fit to the data. 
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Figure 3. The ratio of the 31 GHz excess (extrapolated from the 250 GHz measurement, assuming it is entirely due to free-free emission) 
over the total 31 GHz flux density, against the ratio of the IRAS 60 /im flux density over the 5 GHz flux density. The dashed line is a 
fit to the data that crosses the origin, with a slope of 4.54 it 1.5 and reduced of 0.5. 



6 COMPARISON WITH PREVIOUS WORK 



6.1 NGC 7027 

An important PN which is absent from our sample is north- 
ern NGC 7027, the brightest at radio wavelengths. |Terzian| 
|et al.| ( |1974l ) and |Hoare et al.| ( |1992l ) constructed the SED of 
NGC 7027, finding its radio spectrum follows free-free emis- 



sion. In this r espect NGC 7027 is similar to NGC 6302 (dis- 
cussed in Sec. 8.2.2 1, another high-excitation PN with a mas- 



sive progenitor. Both are surrounded by molecular envelopes 
with copious amounts of dust. Although any cm-excess ap- 
pears to be negligible in NGC 7027 and NGC 6302, the 
mm- wave data could be affected by either molecular lines or 
cold dust, and definitive conclusions in these objects require 
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further analysis. NGC 7027 and NGC 6302 are reminders 
that PNe form a heterogenous class. 



Table 4. Previous flux densities from SEST. Values are reported 
in mjy. Systematic uncertainties are not included, and are as- 
sumed to be of 30%. 



6.2 Previous mm-wave flux densities from the 
literature. 

Deviations from standard free-free spectra in the radio con- 
tinua of PNe were already present in the SEDs obtained 
by |Hoare et al.| ( fl992[ |, although not considered significant. 



PN 



Single pixel 



Heterodyne 



Hoare et al. 



(1992|) that 
show 



It can be seen in Fig. 1 from 
NGC 6572, NGC 6302, NGC 6537 and NGC 6543 
deficits at 200 — 300 GHz. We confirm that our 1.2 mm fiux 
densities are consistent with those measured at 1.1 mm by 
within 2 a, for the four objects we have 



Hoare et al. 



(1992) 
NGC 6572, 



m common: INUU 65Y2, NGC 6302, M2-9, and NGC 6537. 
These objects are studied in detail below. 



6.3 Preliminary mm-wave flux densities from 
SEST. 

The 250 GHz deficit over the expected level of optically thin 
free-free emission was originally found by one of us (L.-A.N) 
with heterodyne and single pixel bolometer observations at 
SEST. But in pointed observations the measurement of flux 
densities is hampered by pointing uncertainties. Until the 
development of bolometer arrays such as SIMBA, designed 
to produce images of the sky, it was difficult to ascertain 
whether the 250 GHz deficits were not due to pointing errors 
(which may also have been the problem in |Hoare et al.|1992[ ). 

The SEST heterodyne observations observations were 
performed between 1995 and 1996 with the 15 m SEST. The 
data were obtained using simultaneously two single-channel 
SIS receivers at 99 and 147 GHz with typical system temper- 
atures of 150 K (SSB) above the atmosphere. The telescope 
beamwidth is 50" at 99 GHz and 34" at 147 GHz (FWHM). 
Two acousto-optical spectrometers (AOS) were used simul- 
taneously with bandwidths of about 1 GHz, channel separa- 
tions of about 0.7 MHz, and resolutions of about 1.4 MHz. 
The telescope was used in a dual beam switch mode with the 
source alternately placed in each of the two beams, a method 
that yields very flat baselines, and allows the determination 
of continuum levels as the offset of the baseline from the 
zero level . The beam separation was about 11^5. The flux 
densities were determined from the chopper- wheel corrected 
antenna temperature, Ta* , using aperture efficiencies, (tia), 
determined from observations of planets, of 0.60 and 0.52 at 
99 and 147 GHz, respectively, and are estimated to be accu- 
rate to within 10 % (one sigma). The pointing was checked 
during the observations using SiO maser sources, and the 
pointing accuracy is estimated to be better than ±3". 

The SEST single pixel bolometer observations were 
performed between 1993 and 1995, using a single channel 
bolometer developed at MPIfR, operating at a wavelength 
of 1200 fj,m (i.e., as for SIMBA). Beam-switching was done 
with a focal-plane chopper having a horizontal beam throw 
of 70". The observations were done in the ON-OFF mode, 
alternating the source between the two beams after 10 s of 
integration time. Pointing wEis checked on nearby quasars 
and was found to be better that 3" rms. The atmospheric 
opacity was determined frequently through skydips. Uranus 
was used as a calibrator source, and the estimated uncer- 
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tainty in intensities are about 10% (one sigma). For single- 
pixel observations we report point-source fiux densities, with 
uncertainties derived from the rms dispersion of all ON-OFF 
pairs. The maps were obtained by scanning the telescope in 
azimuth at a rate of 8 " per second, with adjacent scans 
separated by 8" in elevation. The maps were reduced using 
the NOD2 software ( |Haslam|1974l ) adapted for SEST. Flux 
densities were extracted with Gaussian fits. 

We summarise in Table |4] the measurements from the 
preliminary SEST observations, both from the single pixel 
bolometer and the heterodyne backends. In general the pre- 
vious pointed observations at 230 GHz are within 3 a with 
the 250 GHz values from Table [2] The only exceptions are 
Hb 5, NGC 6537 and NGC 3132, where the pointed data 
gives lower flux densities (consistent with a pointing er- 
ror, or with nebular extensions larger than the beam), and 
NGC 7009, where the difference with the 230 GHz map is at 
3.1 a. The concordance of the whole dataset on NGC 6369 
is particularly comforting, and will be commented further 
below. 



6.4 Absence of radio PN haloes. 

Neither the CBI nor the SIMBA datasets show evidence for 
radio haloes. When compared with the point-source models, 
none of the core-halo models that we tried to fit the CBI vis- 
ibilities gave improved Xcbi values. Extended haloes, about 
4 arcmin in diameter and larger, have indeed been observed 
in a few PNe ( [Corradi et al.||2003[ ). But the optical Ha + 
[Nil] haloes are extremely faint. The halo intensities are of 
order 1/1000 the values in the central regions that dominate 
our radio measurements. The integrated halo line fiuxes are 
about 1/100 the fiuxes of the compact regions. 



6.5 Is the PN cm-wave excess seen in Hll regions? 

As natural comparison objects to PNe, H ll regions are also 
photoionised nebulae with thermal continua. Does their ra- 
dio continua also bear the cm-excess? Hll regions are ex- 
tended and complex objects; the measurement of precise 
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radio spectral energy distributions is difficult. Single-dish 
maps at frequencies higher than ~14 GHz require cancelling 
the sky emission with either fast scanning (as for SIMBA) or 
differencing, while interferometric data are affected by flux 
loss and a wavelength-dependent («, w)-coverage. Another 
difficulty with mm- wave continuum data of Hll regions is 
the possibility of molecular line emission from surrounding 



or coincident molecular material (Sutton et al. 19841 



Notwithstanding these difficulties, [Dickinson et al.| 



( 2007 1 reported that the radio SEDs of six H ll regions con- 



structed from CBl observations and literature data, are com- 
patible with free-free emission. They constrain the possi- 
bility of spinning dust, and note a 3 cr 31 GHz excess in 
G284.3-0.3 similar to spinning dust in its narrow profile. 

But the SEDs of H ll regions at mm- wavelengths are not 
established. A close look at the Boomerang and SEST data 



reported by Coble et al. ( 2003 1 reveals that H ll regions may 



be affected by the same 250 GHz deficit seen in PNe. The 
SED of RCW38 can be constructed with the Boomerang 
data ( |Coble et aL]|2003[ ) and the CBl data poinl[^ We fit 
the 31-400 GHz SED of RCW38 with two power laws. The 
best fit index at high frequencies is a = -1-3.58 ± 0.69, which 
is consistent with grey body emission (with an emissivity 
index of /3 = 1.58 ± 0.69). The low frequency data are fit by 
a = -0.22Io;3g ± 0.14 (where we have indicated the Ax^ = 
1 limits). Given the uncertainties such an index over 31- 
90 GHz is consistent with the indices seen in PNe. 

The H II region SEDs pubUshed by [Malkamaki et al.| 



( 1979 1 also cover the mm-wave regime. The 10-100 GHz 
index is a ~ —0.25 in W51, and is a ~ —0.36 in W3. 
Uncertainties are difficult to estimate from the figures in 



Malkamaki et al. ( |1979) , but our point here is that given 
the existing data, H ll regions may well show the same cm- 
excess as PNe. 



7 CANDIDATE EMISSION MECHANISMS 

In this Section we speculate on the emission mechanisms 
that could explain the observed PN SEDs. We take the case 
of NGC 6369 as an example, since it is the object with the 
most significant 250 GHz deficit. We assume a distance for 
NGC 6369 of 330 pc, and a nebular diameter of 38 arcsec 



(from optical images, as compiled by Acker et al. 1992 1 



Refer to Fig. ^ and Sec. [73] for a comparison of the data 
and a pure free-free spectrum. The free-free continua were 
calculated with the exact Gaunt factors for hydrogen. 

7.1 A synchrotron component? 

Assuming the 250 GHz continuum is due only to free-free 
emission and possibly to the thermal dust Rayleigh- Jeans 
tail, any synchrotron component should be very steep above 
31 GHz, with a spectral index steeper than alf), or about 
-0.4 (Table [l]). This is in contrast with (qI^), which is es- 
sentially flat. 

^ We assigned 10% uncertainties to the CBl data point to allow 
for differences in the photometric extraction procedures: the CBl 
flux densities are obtained by Gaussian fitting, while [Coble et aT] 
l |2003[ l used a circular photometric aperture encompassing most 
of the emission 



Synchrotron self absorption is not expected in PNe. But 
the flat index at low frequencies could be explained by free- 
free absorption with a turn-over frequency at 20-30 GHz. 
For a typical emission measure of 10® cm^pc"^, a turnover 
frequency of 40 GHz requires a free-free screen with an elec- 
tron temperature Te < 600 K for the bulk of the nebula. 



7.1.1 Background 

Synchrotron radiation was sought in PNe during the 1960s 
as a means to detect the magnetic field required by 



Gurzadyan ( 1969 1 in the theory of magnetic shaping (see 
also, e.g., [Pascoli[|T985| |Chevalier fc Luo|[l994[). The evi- 



dence put forward by Gurzadyan ( 1969 1 was based on an 



UV excess in photospheric spectra, and by the mention of 
non-thermal emission being observed in a few objects. How- 
ever, only the case of NGC 3242 appears to be backed by a 



reference. Although Menon & Terzian ( 1965 1 report a non- 



thermal spectrum for NGC 3242, their data is not confirmed 



by the flux densities given in Calabretta ( 1982 I, Condon 



Kaplan ( 1998 1 and Milne & Aller ( 1982 1 



The possibility of non-thermal emission from the inter- 
action of a magnetised and fast wind with dense knots has 
been investigated by Jones et al. ( 1994 1996), in diverse as- 



trophysical contexts. They expect synchrotron emission to 
arise at the shocked interfaces between knots and fast wind. 
Synchrotron should be strongest along the stretched mag- 
netic field lines around the perimeter of the clumps, and 
increase at the time of knot disruption. 

An interesting alternative has been pro posed by[Dgani 



Soker: (1998). In the specific case of PNe, [Dgani fc Soker 



([1998| argue that simultaneous cosmic ray acceleration and 
magnetic field enhancement could occur in the turbulent 
mixing layers at the sides of the knots and their turbulent 
wakes. They find that such synchrotron emission could be 
detectable, at a level of 1 mjy at ~ IGHz, in the case of 
hydrogen deficient nebulae, where the level of free-free emis- 



sion is reduced. Dgani & Soker ( 1998 1 claim the detection in 



PN A 30 of faint (0.1 mJy), compact (3 arcsec), non-thermal 
and circumstellar emission at 8 GHz. But the data have not 
been published. 

There is thus no data available to back the synchrotron 
hypothesis for the cm-excess. Recently [Cohen et ah] ( [2006| 
have reported the detection of two non-thermal hot spots in 
the faint PN surrounding the OH/IR star V1018 Sco (with 
spectral indices of —0.85 ± 0.01 and —0.95 ±0.11 over 5 



and 10 GHz). Cohen et al. (20061 attribute the ionisation of 
this object to the interaction between a fast AGB wind that 
has been recently turned on and is overrunning a precursor 
slow wind. It is very different from the photoionised nebulae 
which make the focus of this work. 

Magnetic fields have however been detected by |Miranda[ 



et al. (2001 \ in the young PN K 3-35, through the detection 



of a circularly polarised OH maser line at 1.665 GHz, imply- 
ing field intensities of order ~ 1 mG. Fields of ~1 mC are 
in the range of the theoretical predictions of [Dgani fc Soker[ 
( 1998 1, based on the equipartition of the magnetic field and 



fast wind energies. 



Vlemmings et al. ( 2006 1 found water masers in the 



proto-PN W43A, whose polarization traces the precession 
of a strong coUimating magnetic fields, > 1 mG on 1000 AU 
scales. Bains et al. (2003 20041 provide further detections 
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of magnetic fields in proto-PNe. Sabin et al. ( 2007 1 observed 



tfie nebular sub-mm dust polarization concomitant to a per- 
vasive ~ f mG field. 

In fact, the energetic requirements of synchrotron radi- 
ation are easily met in PNe. Given a magnetic field and a 
synchrotron index a observed over a frequency range from 
vi to ;/2, the total energy Ue in cosmic rays can be estimated 
from ( |Moffet|p75l ): 



Nj 



2-p 



where, 



Nt = Lobs/ 



3-p 



Et 



3-p 



and 



E 



3eB ' 
1 — 2a, and C3 



1,2, 



2.4 10"'' cgs. 



(1) 
(2) 

(3) 
(4) 



The observed synchrotron luminosity in NGC 6369, between 
0.1 MHz and 31 GHz, is L^bs = 8.7 10^° erg s"\ as obtained 
by integrating a power law spectrum with F^{31 GHz) — 
1.5 Jy and a — —0.3. We find that, for a 1 mG field perme- 
ating the whole nebula, the total energy in cosmic rays re- 
quired to account for Lobs is ~ 5 10*^ erg. Since the lifetime 
of a typical nebula is 10* yr, such an amount of cosmic-ray 
energy requires a stellar luminosity of ~ 4 10""^ L0, much 
less than the observed tip-of-the-AGB luminosity of 10* Lfilf] 
It is interesting to compare the above energetic re- 
quirement and the calculation in Dgani & Soker ( 1998[ ). 
In their idea synchrotron should stem from the whole en- 
semble of globules and their overlapping turbulent wakes, 
encompassing a linear size Rons- Dgani & Soker (19981 as- 



sume equipartition between magnetic field, cosmic ray elec- 
trons, and fast wind. Their Eq. 9 gives a synchrotron lu 
mino sity of only 4.9 1 0^^(-RcnB/5 lO^^^ cm)'^ erg s ^. Accord- 
ingly Dgani & Soker ( 1998p predict synchrotron fiux densi 



ties 1/100 of the values required for the cm-excess. But a 
value of -Rens = 5 10^^ cm is actually very small for a typical 
PN. As seen in the Helix nebula, the ensemble of knots cov- 
ers at least 1/3 of the nebular solid angle, and the nebulae 
considered in this work are typically ~0.1 pc. Given that 



the Dgani & Soker ( 1998 1 synchrotron luminosity scales as 
the cube of Rons, a mere increase to 2.5 10^^ cm, still much 
less than 0.1 pc, approximates the cm-excess. 



7.1.2 NGC 6369 

Given the detection of magnetic fields in PNe, we did the 
exercise of modelling the SED of NGC 6369 including a syn- 
chrotron component (Fig. |4|i). The core free-free emission, 
which represents the bulk of the 250 GHz fiux density, is 
approximated as a uniform slab at 10*K. The synchrotron 
component is modulated by a cold screen, with a power- 
law distribution of opacities (as in the case of M 2-9, see 
Sec. 8.2.1 1, which is required for NGC 6369 to fit the curva- 



ture of the < 30 GHz dat<Q The cold screen opacity includes 
both free-free and bound-free contributions, although the 
latter is negligible at frequencies where the screen opacity is 
~1. 

We could only obtain a fit to the data by including a 
synchrotron break at 31 GHz. For a single power law to 
account for the whole SED, it would need to join the CBI 
and SIMBA data points, so that free-free emission would 
be altogether negligible compared to synchrotron. Thus, the 
best fit we obtained with a single power law has an index of 
—0.57, and a screen at 15K, with a power law distribution 



of opacities (see Sec. 8.2.1 \. This model fits only marginally 
the low frequency data. Note that the need for a high fre- 
quency break is alleviated by the inclusion of a spinning dust 



component (see Sec. 7.2 1 



Any free-free screen with a temperature below ~ 1000 K 
can fit the data, but the electron density in the screen is an 
increasing function of temperature. At a fixed screen tem- 
perature of Tscreon = 500 K, (iVe (scrceu) ) = 1221±291 cm-3 
and iVe(core) = 1446 ± 507 cm"^, while at Tscroon = 10 K 
all of the 250 GHz fiux density is due to the core free- 
free emission, No{coTe) — 2404 cm~'^ and (A'^^ (screen)) — 
33 ± 20 cm~'^. Perhaps an advantage of the lower temper- 
atures is that the low-frequency synchrotron index can be 
somewhat steeper; at 1000 K the index is fixed at ~ —0.1 
by the 5-31 GHz data. 

For integrated flux densities the 2-D opacity fields can 
be modelled following Calabretta ( 1991[ ). Section 



.2.1 



M 2-9 gives more detail on the procedure, as well as an elo- 
quent example of an opacity profile. Fig. [4^ shows the fit 
to the SED of NGC 6369 with a fixed low-frequency syn- 
chrotron index of aLP = —0.4, a fixed high-frequency syn- 
chrotron index of ohf ~ —1.0, a screen at 20 K, wich has an 
average density of 574 ± 88 cm~^ and a power-law opacity 
index 7 — 6.0 (see Sec. 8.2.11. A modified blackbody with 
emissivity index /3 — +2, required to pass through the 100- 
and 60-/im IRAS &ux densities, corresponds to Tdust = 52 K 
and is negligible at < 250 GHz. 



7.1.3 Predictions of the synchrotron interpretation 

The main difficulty for the synchrotron interpretation is thus 
the need for a cold free-free screen. Screen temperatures of 
1000 K are reminiscent of the ORL temperatures, but in 
the bi-abundance model the ORL gas only accounts for a 
very small fraction of the nebular plasma, while the 1000 K 
synchrotron screen represents about half the total nebular 
mass. In contrast, the average electron density of the very 
cold screen is a factor of ~ 1/100—1/1000 less than the bulk 
nebular density. Such a 10-20 K screen could perhaps be 
found in the residual ionisation of dense and largely neutral 
globules exposed to the central star wind, where [Dgani fc] 
Soker ( 1998 1 expect cosmic ray acceleration. 



The cold free-free opacity profile, as a function of a 1- 
D nebular solid angle parameter uj (Sec. 8.2.11, must be a 



^ Synchrotron radiation from tlie nebular knots taps energy from 
tlie fast stellar wind, which in these hot stars is accelerated by 
radiation pressure coupled to the gas through ionic lines (e.g. 
[Pauldrach et alT]|1988[ |. The radiation - fast wind coupling must 
be very effective since observations give P^ad ~ fwind- 



^ This cold free-free component cannot, of its own, explain the 
cm-excess, because the lowest possible free-free index is a = 
—0.175 at 1000 K (see Fig. |A4[ |, and because the nebulae, as 
cold blackbody radiators, would have to be much larger to reach 
the observed < 5 GHz levels 
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power law to attenuate the synchrotron spectrum and mim- 
ick a free-free spectrum in emission below ~31 GHz. Such 
opacity profiles occur naturally in PNe with a large system 
of globules, which seems to be the case in every object ob- 
served at sufficient resolution. A prototypical example is the 



Pacholczyk (1970 his Eq. 3.78). For a synchrotron spectral 



Helix nebula, where knots are known to be dusty ( Meaburn 
et al.|1992 I and molecular ( Huggins et al.|2002[), with phys- 
ical conditions annroachinff the nredictions of'Dvson et al 



ical conditions approaching the predictions ofiPy son et al.| 
(19891: temperatures of T ~ 10 K, proton densities of up 
to ~ 10® cm""^. If the bulk of carbon in the knots is pho- 
toionised by the stellar UV photons with energies < 16 eV, 
then the electron density in globules is about 10^ cm""*. The 
C I free-free continuum could thus provide the cold free-free 
screen required to attenuate in-situ the synchrotron back- 
ground. If so, we except concomitant C I radio recombination 
lines at the sites of synchrotron emission. 

Synchrotron emission should be clumpy as well. The 
cometary globules providing the cold free-free screens must 
be coincident with the synchrotron emitting regions. The 
overlap between cold free-free and synchrotron regions is 
natural in models where synchrotron emission arises in 
the interaction between the fast central star wind and the 



clumps ( [Jones et al. 1996 Dgani & Soker 19981. In the 
clumpy model the free-free opacities should also include a 
filling factor (/, the fraction of nebular volume occupied by 
the clumps). Thus the opacity for a nebula with depth L is 
r = fNeNiUffL, where the free-free absorption coefficient is 
Kg. Values of / ~ 1/10 are allowed by the above synchrotron 
models with a factor of 2-3 increase in the electron density 
Ne- Smaller values of / ~ 1/100 can be obtained with shal- 
lower synchrotron backgrounds, with low frequency indices 
of 0.2. 

If the cm-excess stems from the fast wind, then a rela- 
tionship exists between the cm-excess, the fast wind pres- 
sure Mvoo, and X-ray emission from the shocked cavities 



predicted by the interacting winds model (e.g., Kwok et al 



1978 1 . But the sample of PNe where diffuse X-rays have been 



detected is too small to allow a cross-correlation with our 
sample. Of the PNe listed in Table [2] only NGC 7009 and 
NGC 1360 are among the 13 PNe detected by ROSATPSPC 
( Guerrero et al.||2000 l. XMM and Chandra have allowed to 
resolve the diffuse X-rays and isolate the central star emis- 



sion ( [Guerrero et al.||2005[ ), adding only NGC 3242 to the 
intersection with Table [2[ (an up-to-date list of X-ray PNe 
can be found at http://www.iaa.es/xpn/). 

The data on central star fast winds is also rather scarce, 
and no trends are easily found. For example, NGC 6572, 
with a large cm-excess, has a high mass loss rate of 



M 



y 



in a rather slow wind with Vo 



M = 6 10" 

1190 km s"^ ( [Modighani et al.[[1993[ ). NGC 7009, also with 
a significant excess, has a faster wind 2770 km s ~^ but 



much smaller M = 2.8 10"^ M© yr"^ ( Cerruti-Sola et al 



1989), which gives a ram pressure ~250 times weaker than 



in NGC 6572. IC 418, with a rather small cm-excess, shares 
a similar fast wind as NGC 7009 (vcx, = 940 km s~^ and 
Cerruti-Sola et al.|l98"9 L 



M = 2.8 10"^ M 



yr 



A prediction of the synchrotron interpretation for the 
cm-wave continua of PNe is that the nebular continuum 
should be polarized at ~ 60 GHz. In order to calculate the 
expected level of synchrotron polarization we assume the 
free- free screen is also embedded in a pervasive nebular mag- 
netic field of ImG, and make use of the formulae given in 



index of —0.5, the fraction of intrinsic (unabsorbed) polar- 
ization is n = 69%, in the case of a uniform magnetic field. 
Including a cold nebular screen gives H — 69% at 250 GHz, 
decreasing to H ~ 3 — 15% at CBI frequencies. The total 
nebular emission includes an unpolarized thermal compo- 
nent, which should account for most of the 250 GHz flux 
density to reconcile the continuum and recombination line 
data in the literature. Thus, for a uniform magnetic field, the 
total fraction of polarization Ht should reach a maximum 
of Ht ~ 20% between 50 and 100 GHz. 

But PNe are often very point-symmetric, so that the 
polarization fraction on fiux densities integrated across the 
whole nebulae should vanish for perfect point-symmetry. 



Moreover, in the interpretation of Dgani & Soker ( 1998 1, the 



magnetic fields required for synchrotron emission are ampli- 
fied in turbulent wakes. This is a case of tangled magnetic 
field, without net polarization. Thus, very high angular res- 
olution observations, at 50-100 GHz, are required to detect 
the polarization predicted by the synchrotron interpretation. 

In summary the synchrotron interpretation predicts 
that the radio C I recombination lines should be coinci- 
dent, on angular scales comparable to the cometary glob- 
ules, with peaks in the continuum polarization between 50 
and 100 GHz. 



7.2 Is it spinning dust? 

Spinning dust peaks at ~30 GHz and is 30-70 times weaker 



at 5 GHz (depending on environment, Draine & Lazarian 
1998). If spinning dust emission accounts for 30% of the 



31 GHz flux density, then — 0.045, which is consistent 
with (ai^) within 2 a (see Table [l]). Spinning dust drops 
with a Boltzmann cutoff above 30 GHz, and is > 300 times 
weaker at 100 GHz. Could the 250 GHz deficit be explained 
by a 31 GHz excess due to spinning dust? VSGs are found 
aplenty in PNe, and the proximity of a hot stars warrants 
such VSGs will be charged by the photoelectric effect. 

However, taking all the low-frequency measurements 
into account, we find that the existing models for spin- 
ning dust cannot alone explain the cm-excess in NGC 6369. 
Fig. illustrates that the low-frequency data have too fiat 
a spectral index to be explained by spinning dust. We used 
the spinning dust emissivity of Draine & Lazarian ( 1998 



http : //www . astro . princeton . edu/draine/ ) , for the warm 
neutral medium case (WNM), which peaks at the lowest fre- 
quency. The spinning dust emissivity, being electric-dipolar 
in origin, stems from the Larmor power radiated by individ- 
ual grains, and has a very steep rise with frequency before 
peaking sharply at the Boltzmann cutoff. Magnetic dust. 



proposed by Draine & Lazarian (1999), also with a steep 



low-frequency rise, faces the same problem of not fitting the 
low-frequency data. 

One problem with spinning dust as a candidate emission 
mechanism is that the current models give spectra that are 
too narrow in frequency. Arbitrarily changing the frequency 
peak does not help, since fitting the 5-10 GHz data would 
then miss 31 GHz. The spinning dust emissivities of other 



ISM environments considered by Draine & Lazarian ( 1998 \ 



all peak at slightly higher frequencies. Another problem is 
that the predicted fiux densities, calculated from the highest 
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Figure 4. SED of NGC 6369. We plot flux density in Jy (y-axis) 
as a function of frequency in GHz (x— axis). The CBI and SIMBA 
datapoints are indicated by circles, a- Synchrotron The solid 
line includes a synchrotron component with a break at 31 GHz 
(dash double dotted), a core free-free component at 10* K (dash 
dotted), a sub-mm modified black-body (dotted), and the emis- 
sion expected from the free-free screen, chosen here to be at 20 K 
(dashed), b- Spinning dust The dashed line is a free-free com- 
ponent at 7400 K, with A^e = 2292 cm~^. A spinning dust compo- 
nent for the "WNM" case of |Draine &: Lazarian| | |1998[ l , is shown 
in thick dotted line and contributes to the total flux density shown 
in solid line. The other ISM environments considered by [Draine] 
I&: LazarianI ( |1998^ are shown in dotted line, scaled by the same 
factor as the 'WNM' case. The low frequency data cannot be rec- 
onciled with the spinning dust models, c- Needles The solid line 
includes extinction from very cold metallic needles, the dashed 
line is the free- free component, and the dotted line is the sub- 
millimetre tail of a modified black body, attenuated by gray needle 
extinction. Since needle absorption does not extend appreciably 
below 100 GHz, the dashed line indicates nebular emission had it 
been pure free-free. 



diftuse ISM emissivities (the 'RN' case), barely reach peaks 
of 1 mjy if scaled to nebular densities of 10* cm~'^. 

The PN environment may nonetheless lead to higher 
spinning dust emissivities, especially in the dense cometary 
knots, and the low frequency levels could be explained by a 
synchrotron component. If so, spinning dust, in conjunction 
with a synchrotron component, could explain the cm-excess 
and alleviate the need for a synchrotron break. 



T.3 Extinction due to mm-sized metallic needles? 

Alternatively the 250 GHz deficit could be interpreted as 
extinction due to dust. For standard dust grains such mil- 
limetre extinction would require dust-to-gas mass ratios of 
order 1. However, metallic needles bypass the mass prob- 
lem for extinction, and explain the 250 GHz deficit with 
gray extinction at wavelengths shorter than 1 mm, and an 
opacity decreasing as towards 1 cm. Using the formulae 
given in Dwek ( 2004a|b I , a long- wavelength cutoff for the 



gray extinction of Ao = 1 mm requires a needle aspect ratio 
I /a of ~8000 for a resistivity of pR — 10~® Q cm. Imposing 
a unit needle opacity at 1 mm gives a total needle mass of 
2 10^^ Mq for a grain material density of pm = 8.15 g cm~'^, 
i.e. much less than typical nebular masses of ~ 0.1 M0. Yet, 
to be seen in absorption at 1 mm the needles would have to 
be extremely cold, about 1 K. Thermal balance gives lower 
temperatures for longer needles (Eli Dwek, private commu- 
nication), and although needles are colder than normal in- 
terstellar grains (the needle temperature would be 8 K in 
Gas A, Dwek|2004a \ , temperatures < 1 K seem unphysical. 
A solution may lie in modelling the needle extinction in a 
distribution of needle aspect ratios, such that there exist a 
very small population of ultra-cold grains. 

We attempted to fit the 31-250 GHz drop in NGG 6369 
with needle extinction. The needle extinction is applied as 
a screen to the emergent flux. For simplicity we added the 
needle emission as if it were a foreground to the free-free 
nebula. A fit to the data is obtained with Ao = 1.2 mm, max- 
imum opacity 0.82, and T = 1.2 K. The free-free component 
corresponds to a uniform slab with Ne = 3373 ± 85 cm"'', 
Te = 6580 ± 1031 K. 

The SEST 99 GHz and 147 GHz heterodyne data shown 
on Fig.|4|; are in very good agreement with needle extinction. 
They were originally left out from the fit, and added at a 
later stage for comparison purposes. The SEST heterodyne 
data are subject to numerous uncertainties. The possibility 
of needle extinction will not be considered further in this 
work because it requires extreme temperatures. 



8 COMMENTS ON INDIVIDUAL OBJECTS 

In this Section we fit the SEDs of selected PNe with ana- 
lytic radiative transfer solutions. All fits (except the single 
component free-free spectra) include modified black bod- 
ies constrained to the 60- and 100-/xm IRAS fiux densities, 
with a dust emissivity index of /3 — 2.0. This choice for (3 
may seem rather high compared to the canonical ISM value 
of /3 = 1.7. Yet higher emissivities minimise the contribu- 
tion of the Rayleigh- Jeans tail to the 250 GHz fiux density, 
thereby minimising the 250 GHz deficit relative to the free- 
free level expected from 31 GHz. The PNe listed in Table [2] 
were selected on their being IR-bright. Their flux densities 
peak in the IRAS 25-IRAS 60 bands, and decrease from 60 
to lOO/im. 

Unless otherwise stated the low frequency measure- 

at 5 and 



ments are extracted from Milne & Aller ( 1982 



14 GHz ) |Galabretta| ( |1982| at 408 MHz), and [Condon_& 
|Kaplan| ( [1998| at 1.4 GHz). We generally preferred single- 
dish measurements, when available, in the cases where the 
interferometric data are resolved. PN distances are notori- 
ously difficult to measure, and are usually assigned 100% 
uncertainties. 



8.1 The largest 250 GHz deficits 

8.1.1 NGC 7009 

Few Galactic nebulae can be found near NGC 7009, at 21 h. 
Its CBI and SIMBA flux densities are the best measurements 
reported in this work. The 250 GHz value is averaged from 3 
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the SIMBA data point, at 250 GHz, is low but agrees with 
the Hoare et al. ( 1992 I data within the uncerainties. The 



rising sequence with frequency could perhaps be due to the 
Rayleigh- Jeans tail of sub-mm dust extending to 250 GHz. 
We have also plotted on Fig.[5]data points from the prelimi- 
nary SEST measurements (assigned exagerated error bars). 



8.1.3 NGC 3242 

NGC 3242 has an expansion parallax distance of 420±160 pc 
( [Hajian et al.|1995| , and a nebular diameter of about 25 arc- 
sec. We could fit the SED of NGC 3242 with the same 
synchrotron component as for NGC 6369, except the core 
free-free emission is constrained to the 250 GHz flux den- 
sity, has A^e = 2401 cm~^, and is assigned a temperature of 
20 000 K. A free-free screen at 20 K has an average density 
of 200 cm"'', and a power-law distribution of opacities with 
7 = 10. A single power law fit at 5-31 GHz gives a exponent 



of 



-0.290Io:34o ±0.051, which is different from the free-free 



Qs^ value by more than 3 a. 



Figure 5. SEDs of the PNe with the largest 250 GHz deficits. 
We plot flux density in Jy (y— axis) as a function of frequency in 
GHz (a;— axis). The CBI and SIMBA datapoints are indicated by 
circles. The solid line is a single uniform-slab free-free component, 
drawn to emphasize deviations from free-free spectra, and fit to 
the ^ 31 GHz data with the following parameters: NGC 7009, 
Te = 4550 ± 943 K, = 2618 ± 81 cm'^; NGC 6572, = 
8589±1278 K, Ne = 14810±350 cm"^; NGC 3242, = 7921± 
1601 K, Ne = 2938 ±92 cm'^; NGC 3918, Te = 9750 ±2447 K, 
Ne = 4501 ± 185 cm-3. 



different observing runs. The 31 GHz value is averaged over 
7 different nights and spread over 3 years. 

The SED of NGC 7009 in Fig. |5] is very similar to 
NGC 6369. As in Sec. |7.1| a fit to the data can be ob- 
tained with a broken synchrotron spectrum, with a low fre- 
quency index of qlf = —0.4, and a high frequency index 
a — —1.0 absorbed by a free-free screen with (A^e (screen)) — 
345 ± 54 cm~^, Te (screen) = 50 K, and a power-law opac- 
ity distribution with 7 — 6.0. The core free-free component 
would be Ne = 1456 ± 86 cm'^, Te = 10"* K. We used a 
distance of 420 pc, and a nebular diameter of 28.5 arcsec (as 
compiled by Acker et al.||1992 1 . 



8.1.4 NGC 3918 

The high excitation PN NGC 3918 has been studied in depth 



8.1.2 NGC 6572 

The SED of NGC 6572 in Fig. [5] is again very similar to 
NGC 6369. We adopt the expansion parallax distance of 
1.49 ±0.62 kpc ( Hajian et al.|1995| , and a nebular diameter 
of 8 arcsec (as observed with the VLA by Zijistra et al. 1989 1. 

We could fit the data with an attenuated and broken 
synchrotron spectrum, with a low-frequency spectral index 
of —0.4, and a high frequency index of -1.0. The bulk free- 
free component in Fig. [H] has Ne — 10056 ± 162 cm""' and 
Te = 10* K. We used a free-free screen at 50 K, with 
(A^e (screen)) = 998 ± 248, and power-law opacity distribu- 
tion with index 7 = 5.2. The sub-mm dust has a temperature 
of 76 K. 

The data from [Hoare et al.| ( |1992| | confirm our calibra- 
tion of the SIMBA flux densities. In can be seen in Fig.[5]that 



by Clegg et al. 1 19871. We adopt their preferred distance of 
1.5 kpc, and a nebular angular size of 15 arcsec. The SED 
of NGC 3918 show on Fig. [4] can be fit by a broken syn- 
chrotron spectrum, with spectral indices of —0.4 and —1.0. 
The free- free screen is at 50 K, and has an average density 
of {Ae(screen)) = 371 ± 91 cm'^. The bulk of the 250 GHz 
fiux density is due to the traditional free-free emission, at 
10'' K and Ne = 2342 ± 306 cm'^ 



8.2 Outliers 

In this section we discuss the objects that we have considered 
as outliers from the correlations described in Sec. |5] They 
differ in physical characteristics and radio spectra from the 
rest of the sample. 



8.2.1 M 2-9 

M 2-9 is a peculiar object. Its ionizing central star is under- 
going sustained mass loss, as revealed by broad Ha wings 
(e.g. jArrieta fc Torres-Peimbert||2003| [Swings fc AndrilTat] 



1979fc. As reported by Kwok et al. (19851, the radio SED 



of M 2-9 is that of an optically thick stellar wind, and is in 
marked contrast with the other objects studied here. 

In Fig. [6] we model the radio spectrum of M 2-9 with 
a pow er-law distribution of opacities, following |Calabretta| 
( 1991 1, except the opacity law we used is t{(jj) — 5uj As 
explained in [Calabretta (19911, the detailed spatial distri- 



bution of opacities is irrelevant to the observed free-free flux 
density. Only the distribution of opacities as a function of a 
single solid angle variable is required. In general the opacity 
will be a function of some 1-D parametrisation of solid angle, 
oj. As an example deflnition of a;, |Calabretta, ( ^1991) pixelates 
the sky emission of the nebulae and orders pixels by their 
opacities. A natural choice for uj in spherically symmetric 
nebulae is the nebular radius. In the case of the optically 
thick stellar wind in M 2-9, the bulk of the emission at cm- 
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to mm-wavelengths comes from the wind, with a power law 
profile: t{uj) = 5lj~'' . A clumpy free-free screen can also be 
modelled with a similar opacity profile as M 2-9, except the 
clumps are scattered across the nebula. Larger values of 7 
correspond to smaller fractions of the nebular pixels con- 
taining most of the opacity. Power law profiles allow for a 
smoother transition to the optically thick frequency regime 
than uniform slabs (a proof is given in |Calabretta|1991| ). 
The emergent flux density is 



10 



100 



duj{l — exp{—T^(uj))), 



(5) 



where fiiv is the nebular solid angle, uj = fJ/Ajv is a solid 
angle variable (for instance angular distance from the central 
star), ojti ^ 1 is fixed by requiring dujT{uj) — 1, and 



r,y{uj) = {t„)t{w). 



(6) 



We compiled the SED of M 2-9 with the data in |Hoare et] 



al. ( 1992 1, our own CBI and SEST meeisurements, the 5 GHz 



and 14 GHz data from Milne & Aller ( 1982 1 and the 2.7 GHz 



flux density from Calabretta (19821. We obtain a fit to the 



radio-IR spectrum of M2-9 with a superposition of a uniform 
slab nebula to account for the flattening at low frequencies, 
a power law distribution with S = 1.6 10"'"*, 7 — 2.06 to 
model the optically thick core, an average core density of 
(A^e) = 2093 cm"^, and a temperature Te = 7700 K. The 
best fit average emission measure entering in the average 
opacity (r^) is (EM) = 1.3 lO'^ pc cm~^, with a total nebular 
diameter of 1 arcmin. We take the agreement between the 
SI MBA fl ux d ensity of M 2-9 with th e data from |Hoare et| 
[ah] ( |1992[ ) and [Altenhoff et al.| ( |1994[ ), as well as with the 
model SED, as a confirmation of our data calibration and 
reduction procedure. 



M 2-9 




NGC 6302 




• - - - -J_ 





Figure 6. SEDs of outliers. We plot flux density in Jy (y— axis) 
as a function of frequency in GHz (a;— axis). The CBI and SIMBA 
datapoints are indicated by circles. M 2-9 The components that 
sum up the model spectrum (in solid line) are a modified black 
body (dotted line), an optically thick stellar wind (dashed line) 
and a uniform slab (dash-dotted). NGC 6302 We use a core- 
halo model to fit the observed spectrum. The core component is 
shown in dashed line {Ne = 13470"^, Te = 8917 K), and the halo 
component is shown in dash-dotted line (A^e = 4270 cm~'^, = 
20000 K, with a diameter of 20 arcsec). This model is not unique 
and only serves the purpose of indicating a physical reference. Vy 
2-2 Line styles follow from the NGC 6302 spectrum. The flux 
densities are extracted from [Purton et al.| ( |1982^ and [Altenhoff] 
|et al.||T994l l. 



8.2.2 NGC 6302 

NGC 6302 is the highest excitation PN known, and is prob- 
ably derived from the upper mass range of planetary nebula 
progenitors (e.g., Casassus et al.j2000^ . As such it has been 



the focus of much study. The ionised nebula is surrounded 



by copious amounts of neutral and molecular material ( Ro- 



driguez et al.|1982 1 



In Fig. 
NGC 6302. 



we summarise 
The flux densities 



the available dataset on 
have been extracted from 



Hoare et al. (19921; Rodriguez et aL (19851 



al. (19741; Condon & Kaplan (19981; Zijlstra et al. (19891; 



Terzian et 



Milne & Aller (19821; Milne & Webster (19791; Calabretta 



(19821. We adopt the expansion parallax distance of 1.6 ± 
0.6 kpc (Gomez et al. 19931, and a core nebular diameter 



of 10 arcsec ( Zijlstra et al.|1989 l, surrounded by a halo ex- 
tending up to ~ larcmin. As in the case of M 2-9, we find 
that the SIMBA data point agrees, within the error bars. 



with the data from Hoare et al. ( 1992 1 



The same situation occurs in NGC 6537, often consid- 
ered as a slightly lower excitation twin of NGC 6302. The 
Rayleigh- Jeans tail of the sub-millimetre dust could con- 
tribute to the 250 GHz fiux. 

NGC 6302 is an important cross-check on the SIMBA 
calibration scale. Raising the SIMBA fiux densities by 50% 
would eliminate the cm-excess, and result in a mm-excess 
in some objects. But the SIMBA data on NGC 6302 is con- 
sistent with the literature. Raising the NGC 6302 250 GHz 
fiux densities by 50% could be inconsistent with|Hoare et al.] 



( 1992 1 by 5 <T. 



8.2.3 Vy 2-2 



Vy 2-2 is a very young low excitation PN (e.g., Liu et al 
2001 1, with strong 10/im silicate emission at ~365 K ( Aitken 
fc Rochel|l982) . Faint, compact and opaque radio emission 
at the nucleus of Vy 2-2 is resolved into a 0.4 arcsec diam- 



We fit the SED of NGC 6302 with a core-halo model, as eter ring at 22 GHz (Christianto & Seaquist 19981, and is 



Rodriguez et al. (19851, including the dust temperature 



as a variable. The IRAS 100- and 60-/im data alone give 
Tdust = 46 K, if the dust emissivity is 13 = -1-2. Including 
the UKIRT and JCMT data from iHoare et al.l (|1992|), as 



surrounded by a faint 25 arcsec Ha halo. 

Its distance is greater than 1 kpc (|Christianto & 



Seaquist|1998| ), and traditional indicators place it at 7.9 kpc 
(Zhang 19951 . Since Vy 2-2 is a young object, we adopt 



well as our SIMBA point, gives a dust temperature of 42 K. 
We believe there is room for a colder dust component, that 
could be present at ~200 GHz. 



the distance of 3.9 kpc, given by the IRAS fluxes and the 
tip-of-the-AGB luminosity function ( Casassus et al.|[2001 l. 
Vy 2-2 is very IR bright, while optically thick and faint 
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at radio frequencies. The large uncertainty on its SIMBA 
flux density left it out of Fig. [3] Its IRAS 60 fim to 5 GHz 
ratio is 0.86 10"'', 3-4 times higher than even the largest 
values in Fig. [3] or in NGC 6302. 

The radio SED of Vy 2-2, shown on Fig. |6] can be fit 
with a simple free-free model consisting of a uniform slab, 
with Te = 8522±660, Ne = 3 10^ ±4 10^ cm-^ and a depth 
of 0.01 pc. It can be appreciated this SED is very different 
from the rest of the PNe considered in this work. 

8.2.4 NGC 2440 

NGC 2440 is more evolved than the rest of our sample, 
and its mid-IR flux is due to ionic lines. Dust in this ob- 
ject radiates only in the far-IR, and is presumably located 
further from the central star than in the other PNe. The 
SIMBA data point is unfortunately very uncertain, and so 
is the < 31 GHz data. We mention it here because the low- 
frequency data seems to peak at 5 GHz and drop with fre- 
quency steeper than allowed by free-free emission at low- 
frequencies, in a fashion that is reminiscent of the 250 GHz 
deficit evident in NGC 6369. 



9 CONCLUSIONS 

We have found an inconsistency with the free-free paradigm 
for the radio continua of prototypical PNe, in the form of 
an excess over free- free emission at 31 GHz when compared 
to the 250 GHz fiux densities. The fraction of the 31 GHz 
continuum not due to free-free emission, average over the 10 
objects with less than 30% uncertainties, is 51 ± 3%, with 
a weighted scatter of 11%. The cm-excess is present in all 
objects that are optically thin at 31 GHz, with the exception 
of NGC 6302 and NGC 6537, whose sub-mm dust tails reach 
250 GHz. The most significant detections are in NGC 6369, 
NGC 7009, NGC 3242 and NGC 6572. 

The evidence for the cm-excess stems mostly from the 
250 GHz data acquired with SIMBA. Raising the SIMBA 
calibration scale by ~50% would bring some of the 250 GHz 
levels in agreement with free-free, but would also render the 
data on NGC 6302 and M 2-9 incompatible with the liter- 
ature. Additional evidence can be found in the 5-31 GHz 
spectrum of NGC 3242, and in the average (oii), which are 
steeper than free-free. 

The 31 GHz excess follows a loose correlation with the 
IRAS flux densities. Yet the examination of individual ob- 
jects shows that spinning dust cannot alone explain the low 
frequency data, at least with the models currently available. 
Needle extinction at 250 GHz would require extreme tem- 
peratures. A synchrotron component is a possibility, if it 
is attenuated at low-frequencies by a cold free-free screen 
(< 500 K). The cm- wave - infrared correlation could be due 
to synchrotron radiation originating in compact and dusty 
inclusions exposed to the fast stellar winds observed in PNe. 
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APPENDIX A: THERMAL SPECTRAL 
INDICES 

The free-free emissivity e for electron-ion encounters for a 
dilute plasma (at frequencies much larger than the plasma 
frequency, |Scheuer| |1960') can be written as (e.g. |Beckert"et] 
al.|2000 per unit solid angle). 



8 2n e^Zf 



3 V 3 (mc2)(3/2) VfcT 



exp 



ga, 



(Al) 



where Ni is the ion density and gs{i^, Te,Zi) is a frequency- 
and temperature-dependent Gaunt factor. 
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Al Exact Gaunt factors for free-free emission. 

We searched the hterature for approximations to the exact 
g{f{Zi,i'), and found errors comparable to the level of pre- 
cision in our spectral analysis (of ~5-10%). We thus chose 
to calculate exact Gaunt factors, using the formulae and 



procedures described in Grant ( 1958 1, which have also been 



reproduced in Beckert et al. (20001 and Muiioz (20061. 



The form we used for the non-relativistic Gaunt fac- 
tor in encounters between electrons and bare nuclei is 
that given by [Menzel fc Pekeris] ( |1935[ ) . Velocity-dependent 
Gaunt factors for individual encounters involve the eval- 
uation of complex hypergeometric functions 2Fi{a,b, c, z). 
We used the hypergeometric series when possible (for \z\ < 
min(0.1, |0.1c/(afe)| , as prescribed by Beckert et al.|[2000 |, 
and integrated the hypergeometric equation otherwise. We 
used the GNU Scientific Library routines for ODEs, finding 
that the Bulirsch-Stoer method gave best results. We used 



X 



X^ 




the prescription in Beckert et al. ( 2000 1 as starting point for 
the ODE integration. For \z\ > 1 we used the 'linear trans- 
formation formulae' from Section 15.3 of I Ambramowitz fcl 
Stegun ( 1964 1 to keep the complex argument z within the 



Figure Al. The free-free Gaunt factor for individual electron- 
proton encounters with incident velocities corresponding to 
Tji = —1,-10 and —100. We plot the exact calculation (with 
'x','-|-','*' symbols, respectively), the LFL (dotted line), the clas- 
sical (dashed line), and WKB Gaunt factors (solid line). The re- 
combination limit Gaunt factor (Eq. 6 from [Brussaard Sz van de] 
|Hulst|1962| ) is indicted by a circle. 



unit circle. 

Our hypergeometric routines are robust for any fre- 
quency for relatively fast electrons, when the parameter 
r/i = —aZ{c/vi) < 1, in the notation of [Grant] ( |1958[ ). But 
the ODE integration crashed for large a and b, corresponding 
to slow incident electrons, or when the parameter r/i > 100. 
So we tried the hypergeometric routines in (Press et al.|1996"{ 
NR), edited to double precision and with the starting point 



minus signQ the WKB and the classical Gaunt factor. Oster 



( 1961 1 shows that the LFL converges to the classical limit for 



of Beckert et al. (20001, finding the same numerical limita- 
tion as for GSL. 

The domain of slow incident electrons corresponds to 
the so-called "classical" limit, when the de Broglie wave- 
length is small compared to Z /{mvf), a measure of size 



large rji. It can be seen that as u approaches Vg the exact cal- 
culation scatters about the expected WKB value, especially 
for large rji. This is because the ODE integration becomes 
numerically unstable, and its results cannot be trusted, as u 
approaches Vg. The ODE integration altogether crashes for 
even larger values of rji. But there is always a regime where 
an approximation can be substituted for the exact calcula- 
tion, and this value is the largest between the classical and 
WKB Gaunt factors. Thus, we used the WKB Gaunt factor 
when 1/ > O.Oli^g and rji < —1, and the maximum between 
the classical and WKB Gaunt factors when rji < —100. This 



of the Coulomb field (e.g. Brussaard fc van de Hulst|1962 |. is different from the choice of Beckert et al. (20001, who at 



In what follows we will refer to the "classical" limit as the 
WKB limit, and use the term classical for the classical me- 
chanics calculation leading to the Gaunt factors used by 



Scheuer| (1960) and [Oster] ([1961| . The WKB limit involves 
Hankel functions with imaginary arguments, and are rather 
difficult to evaluate. 

Fortunately, a good and simple approximation for rji > 
100 can be obtained in the following way. When the product 
photon has a frequency close to the recombination (bound- 
free) limit, the WKB Gaunt factor reduces to a very simple 
expression given by Eq. 14 from [Grant[ (]T958). In the limit 
of low-frequencies, the incident electron follows a classical 
trajectory unperturbed by the emission of the radio photon, 
and the Oster/Scheuer Gaunt factors are applicable. At fre- 
quencies between 0.1 and 1000 GHz rji varies between 0.3 
and 2500 when Te = 1 - 50000 K. Fig. |AT] shows the free- 
free Gaunt factor for individual encounters involving inci- 
dent electrons with rji of —1, —10, and —100, as a function 
of frequency. Such values of Vi correspond to the rms veloc- 
ities at temperatures of ~ 10^, 10'^, 10 K, respectively. The 
maximum frequency corresponds to the recombination limit 
fg. We compare the exact calculation, in asterisks, with the 
low- frequency limit (hereafter LFL, as given in, [Beckert et] 
|al.|2000] their Eq. 23, but replacing the absolute value by a ter 1 1961 \ 



ways used the LFL instead of the WKB approximation. The 
relative difference between the use of LFL or WKB is im- 
portant only at very low temperatures, and is within ~ 10"'' 
when Te > 1000 K. 

The Maxwellian integration for the total rate coef- 
ficients (a.k.a. the average coefficients), Eq. 24 in Grant 
( [1958 1, was carried out using the PDL interface to the GNU 
Scientific Library. In order to confirm our calculation we 



(20001 



did the exercise of using the LFL of Beckert et al 
for large \rji\ (i.e. not the WKB limit), and found that we 
obtain the same total Gaunt factors within a relative error 
of ~ 10^^ as the calculations in 



Muiioz 



(20061, who used 



the Mathematica package to directly evaluate thelMenzel & 



20001 



Pekeris ( 1935| ) Gaunt factors (and following |Beckert et al. 



The use of common approximations to the radio- 
frequency Gaunt factors involve errors ~ 5% in the 1- 
300 GHz frequency range for traditional nebular tempera- 
tures. In Fig. |A2|we have compared the exact Gaunt factors 



with the approximations in Scheuer (19601; Oster (19611, 



and the Cloudy photoionisation package (version c06.02c. 



^ The LFL was originally derived by Elwert, as explained by|Os-| 
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Figure A2. Comparison between various gaunt factor formulae 
at Te = 7000 K. We plot the percentual difference between the 
Gaunt factors g for Z = 2 (dotted), the |Oster| | |1961| l (dashed) 
and Cloudy (solid) ( [Ferland et al.|1998| approxi mations, and the 
exact value with Z = 1 from [Menzel &: Pekeris| l |1935[ l , Smp ■ For 
reference we obtain gupif = 250 GHz, Te = 7000 K, Z = 1) = 
2.7954. The curve for Z = 2 has been divided by 2. 




Figure A4. Exact (non-relativistic) optically-thin thermal spec- 
tral indices for encounters between electrons and bare H and He 
nuclei, evaluated for a range of temperatures, and two pairs of 
frequencies: 31-250 GHz, and 5-31 GHz. 




Figure A3. Exact (non-relativistic) Gaunt factors for Z = 1, as 
a function of log(t'/GHz) in x— axis and log(re/K) in j;— axis 



last described by Ferland et al. 19981, which uses free- free 
Gaunt factors based on Hummer (19881. Note that the Al- 



tenhoff et al. ( 1960[) free-free opacities, with the a{u,T) fac- 



tor correction from |Mezger fc Henderson ( 1967 1, are equal to 
free-free opacities calculated directly from the Oster ( 1961 1 
Gaunt factors. 

For practical applications, in order to speed up the op- 
timisations involved in fitting the PN SEDs, we gridded the 
exact free-free Gaunt factors in the (log(re), log(i')) space. 
Fig. |A3| shows the corresponding Gaunt factor images. It can 
be appreciated that there are very little differences between 
the Z — 2 and Z = 1 Gaunt factors. The limit for a cold 
plasma, at high frequencies, is 1, as expected. Note that the 
Gaunt factors given by Oster ( 1970| tend to for cold plas- 
mas because they include the exp(— /ii//(fcr)) term, which is 



usually associated with the emissivity (Eq. Al I rather than 
with the Gaunt factor. 



into accounlpl W e use the formulae compiled in Brussaard 
& van de Hulst (19621, taking free-bound Gaunt factors of 
1. Fig. 8 of Brussaard & van de Hulst ( 1962 I shows that, in 



the radio-sub-mm regime, the free-bound Gaunt factors are 
essentially 1 for recombinations to any principal quantum 
numbers. 



A3 Contribution from He and representative PN 
spectral indices. 



The free-free a nd f ree-bound emissivities depend on temper- 
ature. In Figs. 

and aQi"(He^"'"), in the case of 



A4 w 

(He++), a!f(H 
optically-thin emission and for a range of electron tempera- 
tures Te. 

The Gaunt factors for the He free-free continuum should 
be almost identical to the H values in the domain where 
the wavelength of the electron incident on He"*" has a wave- 
length much larger than the He^ Bohr radius. Thus, we 
can neglect differences with He up to a temperature Te = 
(2/i/ao)^/(3fcme) « 17 10*^ K. As explained above, the 
Gaunt factors for the He"*" free-free continuum are slightly 
different from H. Considering that the most He rich nebulae 
(of Peimbert type I) reach He abundance by number of 18% 
relative to H, a fraction of which is doubly ionised, we can 
conservatively assume a contribution from the He^ contin- 
uum of up to 20%, and a minimum of zero. This translates 
into the spectral index uncertainties given in Table [l] where 
we have assumed an average abundance of He^"*" of 10% 
relative to H"*". 



A2 Free-bound emission 

^ we have not included the contribution of radio recombination 
Since we are considering thermal emission from cold plasma, hues (RRLs) to broad band data. RRLs could be important at 

the possibility of radio free-bound emission must be taken ~300 GHz for Te --^ 20 K 



